/ 

j 



DOUBLY CONSTRAINED ROBUST CAPON BEAMFORMER 



Inventor(s): Jian Li 

Petre Stoica 
Zhisong Wang 



UNIVERSITY OF FLORIDA 



(UF Disclosure No. 11190) 



AKERMAN SENTERFITT DOCKET No.: 5853-417 



Express Mail No. EV346756443US 



{WP148440;2} 



DOUBLY CONSTRAINED ROBUST CAPON BEAMFORMER 



CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] Not applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

[0002] The United States Government has rights in this invention pursuant to National 

Science Foundation grants CCR-0104887 and ECS-0097636. 

FIELD OF THE INVENTION 
[0003] The invention relates to Capon beamforming, and more specifically to an advanced 
Capon beamforming method and related apparatus which provides computational efficiency and 
power estimates for signals of interest which can be adjusted for steering vector uncertainty. 

BACKGROUND 

[0004] Beamforming is a ubiquitous task in array signal processing, such as for radar, sonar, 
acoustics, astronomy, seismology, communications, and medical imaging. The standard data- 
independent beamformers include the delay-and-sum approach as well as methods based on 
various data-independent weight vectors for sidelobe control. The data-dependent Capon 
beamformer adaptively selects the respective weight vectors to minimize the array output power 
subject to the linear constraint that the signal of interest (SOI) does not suffer from any distortion 
(unity gain and no phase shift). 
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[0005] The Capon beamformer has better resolution and better interference rejection 
capability as compared to the data-independent beamformer, provided that the array steering 
vector corresponding to the SOI is accurately known. However, in practice the SOI steering 
vector is assumed, rather than being accurately known. This results in steering vector error. 
Steering vector error generally occurs because of differences between the assumed SOI arrival 
angle and the true SOI arrival angle and between the assumed array response and the true array 
response, such as due to array calibration errors. Whenever the SOI steering vector error 
becomes significant, the performance of the Capon beamformer can become worse than standard 
beamformers. 

[0006] Many approaches have been proposed during the past three decades to improve the 
robustness of the Capon beamformer. To account for array steering vector errors, additional 
linear constraints, including point and derivative constraints, can be imposed. However, these 
constraints are not explicitly related to the uncertainty of the array steering vector. Moreover, for 
every additional linear constraint imposed, the beamformer loses one degree of freedom (DOF) 
for interference suppression. It has been shown that these constraints belong to the class of 
covariance matrix tapering approaches. 

[0007] Diagonal loading (including its extended versions) has been a popular approach to 
improve the robustness of the Capon beamformer. The diagonal loading approaches are derived 
by imposing an additional quadratic constraint either on the Euclidean norm of the weight vector 
itself, or on its difference from a desired weight vector. Sometimes diagonal loading is also 
proposed to alleviate various problems of using the array sample co-variance matrix and to better 
control the peak sidelobe responses. However, for most of these methods, it is not clear how to 
choose the diagonal loading based on the uncertainty of the array steering vector. 
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[0008) The subspace based adaptive beamforming methods require the knowledge of the 
noise covariance matrix. Hence they are sensitive to the imprecise knowledge of the noise 
co variance matrix in addition to the array steering vector error. Making these methods robust 
against the array steering vector error will not cure their sensitivity to imprecise knowledge of 
the noise covariance matrix. 

[0009] Most of the early suggested modified Capon beamforming algorithms are rather ad 
hoc in that the choice of parameters is not directly related to the uncertainty of the steering 
vector. Only recently have some Capon beamforming methods with a clear theoretical 
background been proposed, which, unlike the early methods, make explicit use of an uncertainty 
set of the array steering vector. However, even Capon methods which make explicit use of an 
uncertainty set of the array steering vector are computationally inefficient and cannot generally 
provide accurate power estimates for signals of interest. 

[0010] An improved Capon beamformer method referred to as a robust Capon beamformer 
(RCB) is described in copending and commonly assigned U.S. Application No. 10/358,597 
entitled "Robust Capon Beamforming" by the same inventors as the current application. U.S. 
Application No. 10/358,597 is hereby incorporated by reference into the current application in its 
entirety. The method described therein includes the steps of providing a sensor array including a 
plurality of sensor elements, wherein an array steering vector corresponding to a signal of 
interest (SOI) is unknown. The array steering vector is represented by an ellipsoidal uncertainty 
set. A covariance fitting relation for the array steering vector is bounded with the uncertainty 
ellipsoid. The matrix fitting relation is solved to provide an estimate of the array steering vector. 
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SUMMARY OF THE INVENTION 

[0011] A method for enhanced Capon beamforming is referred to herein as "doubly 

constrained robust Capon beamforming ,, (DCRCB). The method includes the steps of providing 
a sensor array including a plurality of sensor elements where an array steering vector 
corresponding to a signal of interest (SOI) is unknown. The sensor elements can comprise 
antennas or ultrasound transducers. 

[0012] A covariance matrix fitting relation is bounded for the unknown array steering 

vector by constraining the array steering vector using a constant norm constraint and a spherical 
uncertainty set constraint. Although a more general ellipsoidal uncertainty set can theoretically 
be used with the invention, the problem becomes substantially more difficult to solve as 
compared to a spherical uncertainty set. A matrix fitting relation is then solved to provide an 
estimate for the array steering vector. The Lagrange Multiplier Method is preferably used to 
solve the matrix fitting relation. 

[0013] The method can include the step of determining the power of the SOI, including 

adjustment for the estimated array steering vector. The SOI power can then be used to determine 
a direction of arrival (DOA) of the SOL The method can also be used to determine a weight 
vector for said sensor array. The method can include the step of determining a waveform of the 
SOI from the weight vector. 

[0014] A sensor-based system includes a sensor array including a plurality of sensor 

elements, wherein an array steering vector corresponding to a signal of interest (SOI) is 
unknown. The system includes a signal processor communicably connected to the elements of 
the sensor array. The signal processor bounds a covariance matrix fitting relation for the array 
steering vector by constraining the array steering vector using a constant norm and a spherical 
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uncertainty set, and solving said matrix fitting relation to provide an estimate of the array 
steering vector. The sensor system can be a radar, cellular communications, sonar or acoustic 
imaging system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] A fuller understanding of the present invention and the features and benefits thereof 
will be accomplished upon review of the following detailed description together with the 
accompanying drawings, in which: 

[0016] FIG. 1 is a diagram of an ultrasound imaging system which includes an advanced 
signal processor, according to an exemplary embodiment of the invention. 

[0017] FIG. 2(a) and (b) are power estimates and diagonal loading levels (using R) versus 
the steering direction 6 when e = 1.0 and j3 = 6.0. The true powers of the incident signals from 

-35°, -15°, 0°, 10°, and 40° are denoted by circles, and e 0 = 1.0. 

[0018] FIG. 3(a) shows for each steering direction 0, j3 being chosen to make NCCB have 
the same diagonal loading levels as RCB with e = 1 .0. 

[0019] FIG. 3(b) are power estimates versus the steering direction 6 via the RCB and 
NCCB approaches. For the RCB, e = 1 .0. For the NCCB, |8 is chosen as in FIG. 3(a). The true 
powers of the incident signals from -35°, -15°, 0°, 10°, and 40° are denoted by circles, and 

€o-1.0. 

[0020] FIG. 4(a) and (b) are power estimates versus the steering direction 6 when (a) e = 0.7 
and (b) e = 1.5, respectively. The true powers of the incident signals from -35°, -15°, 0°, 10°, 
and 40° are denoted by circles, and e 0 =1.0. 

[0021] FIG. 5(a) - (j) compares the RCB and DCRCB approaches for each incident signal, 
as e varies, when e 0 = 1 .0. 

[0022] FIG. 6(a) - (d) compares the RCB and DCRCB approaches with (a), (b) e = 1 .0 and 
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(c), (d) e = 2.0 when 0 4 (the direction of arrival of the fourth signal) is changing from 20° to 60° 
The SOI power is 10 dB and e 0 = 1 .0. 

[0023] FIG. 7(a) - (d) compares the RCB and DCRCB approaches, as the snapshot number 
varies, when (a), (b) e = 0.6 and (c), (d) e = 2.0. The SOI power is 10 dB and e 0 = 0.5603 
(corresponding to A = 2.0°). 

[0024] FIG. 8(a) - (£) compares the RCB and DCRCB approaches with (a), (b), (e) e = 0.6 
and (c), (d), (f) e - 2.0. The SOI is 10 dB and e 0 = 0.5603 (corresponding to A = 2.0°). 
[0025] FIG. 9(a) - (f) compares the RCB and DCRCB approaches as SNR varies when (a), 
(b) e 0 = 0.1474, (c), (d) e 0 = 0.5939, and (e), (f) e 0 = 1.2289, corresponding to A = -1.0°, -2.0°, 
and 3.0°, respectively. The SOI power is 10 dB and e = 0.6. 

[0026] FIG. 10(a) - (f) compares the RCB and DCRCB approaches in 100 Monte-Carlo 
trials when (a), (b) e = 0.6, (c), (d) e = 2.0, and (e), (f) e = 3.0. The direction of arrival of the 
fifth signal is uniformly distributed between 38° and 42° and its assumed angle is 40°. The SOI 
power is 10 dB and 0 < e 0 < 0.5939. 

[0027] FIG. 1 1 (a) - (d) compares the RCB and DCRCB approaches in the presence of 
potentially significant array calibration errors. The SOI power is 10 dB, and 0 < e 0 < oo . 
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DETAILED DESCRIPTION 

[0028] The invention relates to an improved beamforming method and related apparatus. A 
method of doubly constrained robust Capon beamforming (DCRCB) includes the steps of 
providing a sensor array including a plurality of sensor elements, where an array steering vector 
corresponding to a signal of interest (SOI) is unknown. As used herein, an "unknown" steering 
vector refers to a condition where the array steering vector is at least not precisely known, and 
may not even be accurately known. 

[0029] A covariance matrix fitting relation is bounded for the unknown array steering vector 
by constraining the array steering vector using a constant norm constraint and a spherical 
uncertainty set constraint. The matrix fitting relation is then solved to provide an estimate of the 
array steering vector. The DCRCB is a natural extension of the standard Capon beamformer 
(SCB), obtained via covariance matrix fitting, to the case of uncertain steering vectors by 
enforcing the constraints described above. 

[0030] The RCB described in U.S. Application No. 10/358,597 and the DCRCB described 
herein are computationally more efficient and conceptually simpler than other RCBs. Improved 
computational efficiency which can increase processing speed and reduce hardware cost over 
existing methods and beamforming based systems. 

[0031] Both the RCB described in U.S. Application No. 10/358,597 and the DCRCB 
described herein address the scaling ambiguity problem, which can be very important for SOI 
power estimation, while the other RCBs ignore it. The difference between DCRCB described 
herein and the related RCB disclosed in copending U.S. Application No. 10/358,597 entitled 
"Robust Capon Beamforming" by the same inventors as the current application lies in one 
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additional constant norm constraint on the array steering vector in the former. This constraint 
results in the quite distinct performances between DCRCB and RCB. Examples provided herein 
demonstrate performance differences between the DCRCB and the RCB under several selected 
conditions. Generally, for a reasonably tight spherical uncertainty set of the array steering vector, 
the DCRCB described herein is preferred to the RCB in U.S. Application No. 10/358,597 for 
applications requiring high SINR, while the RCB in U.S. Application No. 10/358,597 is 
generally favored for applications requiring accurate signal power estimation. 

[0032] Although this specification primarily addresses estimation of the array steering 
vector (a 0 ) for the SOI, the DCRCB is also described for SOI waveform estimation, direction of 
arrival (DOA) estimation, source localization and imaging, SOI power estimation, and weight 
vector determination. 

[0033] • To derive the DCRCB, a reformulation of the Capon beamforming problem is used, 
to which the uncertainty set ( || a 0 - a || 2 <C) and the norm constraint ( | a | 2 = M) is appended. 
Proceeding in this way a robust estimate of the signal of interest (SOI) power is directly 

obtained, without any intermediate calculation of a weight vector w . 

max a 2 subject to R - a 2 aa* > 0 

<r 2 ,a 

Ha - 5|| 2 < e 
INI 2 = M 
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The DCRCB problem can be reduced to the following problem: 

mina*R _1 a subject to ||a - a|| 2 < e 

a 

||a|| 2 = M 

using l|a|| 2 = ||a|| 2 = M 

the following results: 

mWR^a subject to Re(a*a) > M - \ 

||a|| 2 = M 

Consider first the problem without the uncertainty set: 

mina^R^a subject to ||a|| 2 = M 

a 

The solution to the above problem is the eigenvector of R corresponding to its largest 
eigenvalue, scaled so that 

||a|| 2 = M 

However, the eigenvector of a matrix is unique only up to a scalar. Even with the normalization 
above, the phase of must be chosen so that &■ &■ is real-valued and non-negative and 
hence R e O*a) is maximum . if it satisfies 116 - M ~ 2, then it is the desired 

solution to the DCRCB problem and the uncertainty set is an inactive constraint. 
If not, i.e., 
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Re(a*a) < M - c/2 

and the following problem is solved: 

/i(a, A, A) = a*R" x a + A (||a|| 2 - m) + /t (2M - e - a*a - a*a) 

which can be written as: 

MaXfi) = [«-/i{R- 1 + ^- 1 *]*(R- l + ^[a-/i(R- l + Xl)- 1 a] 

-£ 2 a* (R- 1 + Al) -1 * - AM + /2(2M - e) 
The unconstrained minimizer for the above function is given by: 

a^ = MR -1 + Air's 

Thus: 

/ 2 (a^,A,/2) = -// 2 a*(R _1 +AI) _1 a-AM + /i(2W-e) 
< a*R _1 a for any a € <S 

Differentiation of the above function with respect to ^ gives 

2M-e 
^ = 2a*(R- 1 + AI)" 1 a 

(M-f) 2 

/3(^,A,A) = -AM + _^_ 1+ ^_ la 



Differentiation of the above function with respect to ^ gives: 
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h(X) = p 



where: 



a^R" 1 +AI)- 2 a 
a^R-i + Al)- 1 ^] 



h{\) = - 



and: 



p= 



M 



It can be shown that M}) is a monotonically decreasing function of A , which can be obtained 
efficiently via several methods, such as Newton's method. Then, the following is obtained: 



e\ (R^ + AI)- 1 * 



a*(R~ a + AI)- 1 ^ 



and: 



w = 



a*R-!a 



f / i \-i 1 (r+ti) a 

V 2J I \ A/ Ja*(R+il) R(R+il) a 



-l 



[0034) The DCRCB weight vector (w ) has the form of diagonal loading except for the real- 
valued scaling factors as well as the diagonal loading level being possibly negative. However, 
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the scaling factor is not generally important since the quality of the SOI waveform estimate is 
typically expressed by the signal-to-interference-plus-noise ratio (SINR). 



SINR = 0 1 - - 1 - 



w* (Eti^^ + Q)w 

[0035] Due to the norm constraint on the steering vector, the constraint set in the DCRCB 
problem is no longer convex. To focus on the convex optimization problem, which is considered 
to be one of the simplest classes of optimization problems, the RCB algorithm was developed. 
RCB is an approximate solution to the DCRCB problem by first finding a solution for the 
problem by ignoring the norm constraint and later imposing it on the solution. More specifically, 
the RCB is derived by first solving the following optimization problem: 

mina*R~ 1 a subject to ||a - a|| 2 < e 

a 

After a is obtained, the norm constraint 1 1** 1 1 —M 
is imposed when estimating the SOI power as: 

,2 ||a" 2 
a 0 = 



Ma*R _1 & 

a*U(I + Xr)~ 2 r 2 u*a 
Ma*u(i + Ar)- 2 ru*a 



[0036] Even though the RCB is an approximate solution to the original problem, numerical 
examples presented in Application No. 10/358,597 demonstrate its excellent performance. 
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Moreover, the spherical uncertainty can be readily generalized to both non-degenerate and flat 
ellipsoidal uncertainty sets for RCB. It is, however, not as easy to generalize to the case of 
ellipsoidal uncertainty sets for DCRCB, which requires a two-dimensional search to determin. 
the Lagrange multipliers. 

[0037] In many applications, such as in communications or global positioning, the focus 
on SOI waveform estimation. As demonstrated below, the weight vector can first be obtained 
from the steering vector estimate provided by the DCRCB. The weight vector can then be use 
to weigh the array output to obtain the desired SOI waveform estimate. 

Let s 0 (n) denote the waveform of the SOI. Then, once the estimated SOI steering ve 
is obtained with the DCRCB, s 0 («) can be estimated as in the SCB as follows: 



So(n) = w * Xl 



n 



where a in: 




is used in lieu of a 0 in 



w 0 = K a» 



a 0 * R" 1 a G 



to obtain W : 



W = 



R -1 a 
a*R -1 a 



( R+ i I )" la 



-1 



Ra 
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Note that the DCRB weight vector ( W) also has the form associated with the diagonal loading 
approach, except for the real-valued scaling factor in the equation above as well as the fact that 
the diagonal loading level in this equation can be negative. However, the scaling factor is not 
generally significant since the quality of the SOI waveform estimate is typically measured by the 
signal-to-interference-plus-noise ratio (SINR): 

SINR = fo\***>\ 2 

w*(£f=i*Ka£ + Q)w 

which is independent of the scaling of the weight vector ( W). 

[0038] The DCRCB can also be used to obtain the power estimate corresponding to each 
direction angle or each coordinate in a two- or three-dimensional space. The power estimates can 
be used for DOA estimation, source localization, and imaging. The power estimates as a 
function of the two- or three-dimensional coordinate are two- or three-dimensional images, 
respectively. 

[0039] A wide variety of improved systems can be based on the invention, particularly 
wherever high signal-to-interference-plus-noise-ratios (SDSfRs) and/or accurate power estimates 
are desired and there is some uncertainty in the array steering vector and/or the snapshot number 
is limited. The knowledge of the SOI steering vector can be imprecise, which is often the case in 
practice. Whenever this happens, the performance of the standard Capon beamformer may 
become worse than that of the standard data-independent beamformers. 

[0040] The DCRCB can be used to cure this problem and the power estimates are much 
more accurate. In addition, high signal-to-interference-plus-noise-ratios (SINRs) can be 
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achieved. Furthermore, the DCRCB can work even when the number of snapshots is less than 
the number of the elements in an array. Thus, the invention is expected to directly replace some 
data-independent beamforming methods which have marginal performance for interference 
suppression, and some data-dependent methods which are sensitive to the steering vector 
mismatch. 

[0041] In the exemplary embodiment shown in FIG. 1, a 3D ultrasound imaging system 100 
according to the invention is shown. An object 20 is irradiated by a source 64 which emits 
ultrasound. Structures 30 within the object 20 reflect the ultrasound waves in different directions. 
The reflected ultrasound waves are captured by a plurality of discrete sensors 80 of 2D or 3D- 
sensor array 60. Output signals of the sensors 80 are transmitted via a suitable communication 
link 40 to a processing unit 42, such as a conventional computer workstation for processing. 
Reconstructed images can then be displayed by display 44. High resolution images and compact 
design of the multidimensional sensor array 60 and the processing unit 42 are made possible 
using the advanced DCRCB process according to the present invention. 

[0042] Improved adaptive array antennas based on the invention can be used in various 
mobile communication systems including base-mobile, indoor-mobile, satellite-mobile, and 
satellite-satellite communication systems. They can help improve the system performance by 
increasing channel capacity and spectrum efficiency, extending range coverage, tailoring beam 
shape and steering multiple beams to track a large plurality of mobile units. Furthermore, such 
systems can reduce multi-path fading, co-channel interference, system complexity and cost, as 
well as the bit error rate. 

[0043] The invention can also be used for improved Global Positioning Systems (GPS). In 
GPS applications, the navigation signals transmitted from satellites generally have low power 
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and very strong interferences may co-exist with the navigation signals. In addition, multi-path is 
introduced by the environment. Consequently, severe performance degradation may occur. 
Adaptive arrays according to the invention can be used to achieve a high signal-to-interference- 
plus-noise ratio (SINR) and accordingly improve the positioning accuracy of GPS. 

[0044] The invention can also be used for improved Quadrupole Resonance (QR) explosive 
detection and Magnetic Resonance Imaging (MRI). Improved adaptive arrays can also be used 
to obtain high-resolution images for the target detection and parameter estimation as well as 
interference suppression. Adaptive arrays can also be used to obtain more accurate power 
estimates of acoustic sources in the presence of interferences and noise within a reverberant 
environment. 

[0045] The invention can also be used for improved radar and sonar. Regarding radar and 
sonar, very strong jammers and clutter may be present. Adaptive arrays according to the 
invention can be used to improve the SINR and obtain high-resolution imaging results. 

[0046] The invention can also be used for improved speech enhancement. For example, in a 
conference room, the desired speech signal from the mouth of a talker may be corrupted by 
interfering signals such as other talkers and room reverberation. By exploiting adaptive arrays 
based on the invention, a higher-quality speech signal can be obtained. 

[0047] The invention can also be used for speech source localization and tracking. In this 
embodiment, improved adaptive arrays can be used to locate and track one or more speaking 
sources. The speaker location estimates can be employed to steer a camera or a series of cameras 
and frame the active speaker in a video-conference system. 
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[0048] The invention can be used for improved hearing aids. Adaptive arrays can be used to 
reduce the effect of background noise and interferences to significantly improve the speech 
understanding of the hearing impaired and to increase their overall satisfaction with the hearing 
aid device. 

[00491 The invention also has applications with automobiles. For example, in many 
countries and regions hand-held telephony in cars is prohibited by legislation. The advantages of 
hand-free telephones in cars are safety and convenience. However, by installing the microphone 
far away from the user, poor sound quality and acoustic echoes may occur. Improved adaptive 
arrays according to the invention can be used in cars for improved speech enhancement, echo 
cancellation and reverberation suppression. 



Examples 

[0050] It should be understood that the examples and embodiments described herein are for 
illustrative purposes only and that various modifications or changes in light thereof will be 
suggested to persons skilled in the art and are to be included within the spirit and purview of this 
application. The invention can take other specific forms without departing from the spirit or 
essential attributes thereof. 

[0051] Numerical examples are provided in this section to compare the performances of 
the delay-and-sum beamformer, SCB, NCCB, RCB and the DCRCB according to the present 
invention. In all of the examples considered below, a uniform linear array with M= 10 sensors 
and half-wavelength sensor spacing is assumed having a spatially white Gaussian noise 
co variance matrix given by Q = I. 
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[0052] For NCCB, = M' J/r ^ aeK P (P ^ -0 is set as a user parameter. The larger 
the 0, the closer NCCB is to SCB. On the other hand, the smaller the ft the closer NCCB is to the 
delay-and-sum beamformer. When 0 = 1, the NCCB becomes the delay-and-sum beamformer 
and hence it uses the assumed array steering vector divided by M as the weight vector. The true 
steering vector error s 0 = min a \\z 0 e Ja - if is defined, where a 0 is the true steering vector and 
a is the assumed one. Unless otherwise stated, the beamforming methods are used with the 
theoretical array covariance matrix R. 

[0053] First, an imaging example is considered where the incident signal power as a 
function of the signal arrival angle 8 relative to the array normal is to be determined. It is 
assumed that there are five incident signals with powers 30, 60, 40, 35, and 10 dB from 
directions -35°, -15°, 0°, 10°, and 40°, respectively. To simulate the array calibration error (the 
sensor amplitude and phase error as well as the sensor position error), each element of the 
steering vector for each incident signal is perturbed with a zero-mean circularly symmetric 
complex Gaussian random variable normalized so that e 0 =1.0. The perturbing Gaussian random 
variables are independent of each other. For RCB and DCRCB, e = 1.0 is used. For NCCB, 0 = 
6.0 is chosen so that the peak widths of the NCCB and DCRCB are about the same. Figure 2(a) 
shows the signal power estimates as functions of the arrival angle 6 obtained via using the delay- 
and-sum beamformer, SCB, NCCB and DCRCB methods. The small circles in the figure denote 
the true (direction of arrival, power)-coordinates of the five incident signals. Since the power 
estimates of RCB and DCRCB are almost the same for this example, only the DCRCB power 
estimates are shown in the figure. Note that SCB can give good direction-of-arrival estimates for 
the incident signals based on the peak locations. However, the SCB estimates of the incident 
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signal powers are significantly off. The NCCB is more robust than SCB but still substantially 
underestimates the signal powers. On the other hand, the DCRCB according to the invention 
provides excellent power estimates of the incident sources. As expected, the delay-and-sum 
beamformer has poorer resolution than the other beamformers. Moreover, the sidelobes of the 
former result in false peaks. 

[0054] Figure 2(b) shows the diagonal loading levels of the NCCB, RCB and DCRCB 
approaches. The NCCB can have a nonzero or zero diagonal loading level. This results in the 
discontinuities in the NCCB diagonal loading level curve. The discontinuity in the DCRCB 
diagonal loading level curve is shown always around the strongest signal. As a result, DCRCB is 
no longer a diagonal loading approach around the strongest signal. 

[0055] For each steering angle 6 in Figure 3(a), /3 is chosen to make NCCB have the same 
diagonal loading level as RCB when e = 1.0 is used in RCB. Note that for NCCB and RCB to 
have the same diagonal level, 0 must be chosen in a complicated manner depending on both e 
and the data itself. Figure 3(b) shows the signal power estimates as functions of 6 obtained via 
using NCCB and RCB with the 0 in NCCB chosen so that NCCB and RCB have the same 
diagonal loading levels. It is noted that the RCB signal power estimates are much more accurate 
than those obtained using NCCB and hence the norm constraint imposed on a is very helpful for 
accurate SOI power estimation. 

[0056] Figures 4(a) and 4(b) show the power estimates as functions of 6 obtained via using 
RCB and DCRCB with e = 0.7 and e = 1.5, respectively, for the example in Figure 2. Note that 
when e < e 0 = 1 .0, the RCB and DCRCB signal power estimates are not as accurate as when 

c > e 0 , but the peaks are sharper. 
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[0057] In Figure 5, the SDSfRs and the signal power estimates for the five incident signals 
are compared, as functions of e ? obtained via using RCB and DCRCB. Figures 5(a), 5(c), 5(e), 
5(g), and 5(i) show the SINRs of the five signals as functions of e. Figures 5(b), 5(d), 5(f), 5(h), 
and 5(j) show the power estimates of the five signals as functions of e, with the horizontal dotted 
lines denoting the true signal powers. Note that except for the 4th signal, the SINR of DCRCB is 
in general higher than that of RCB when e is not too far from e 0 . Hence for applications requiring 
waveform estimation, the former may be preferred over the latter if Go is known reasonably 
accurately. For the 2nd signal in Figure 5(c), the DCRCB uses the scaled principal eigenvector as 
the estimated steering vector. For this case, DCRCB is always better than RCB, no matter how e 
is chosen. On the other hand, for signal power estimation RCB in general outperforms DCRCB 
and hence may be preferred in applications such as acoustic imaging where only the signal 
power distribution as a function of angle or location is of interest. Also the larger the e, the more 
RCB and DCRCB will overestimate the signal power. Therefore, if possible, e should not be 
chosen much larger than eo. 

[0058] In the next examples, the 5th signal from 40° is treated as the signal of interest (SOI). 
The other 4 signals are considered as interferences. In the following figures for the SOI power 
estimates, the dotted lines correspond to the true SOI power. 

[0059] First, a scenario is considered where the fourth signal changes its direction of arrival 
from 20° to 60° with the directions of arrival of the SOI and the other three interfering signals 
fixed. The array suffers from the same calibration error as in Figure 2. Note from Figure 6 that 
when the direction of arrival of an interference signal becomes too close to that of the SOI, both 
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RGB and DCRCB suffer from severe performance degradations in both SINR and SOI power 
estimation accuracy. As expected, the larger the e used the weaker the interference suppression 
capability of both methods when an interfering signal is nearby the SOI. 



[0060] The effect of the number of snapshots N on the SINR and SOI power estimation 
accuracy of RCB and DCRCB are considered when the sample covariance matrix R, from 

1 N 

n=,l 

is used in lieu of the theoretical array covariance matrix R. It is assumed that the steering vector 
error is due to an error in the SOI pointing angle, which is assumed to be 0 5 +A, where 8 5 is the 
true arrival angle of the SOI. In this example, c 0 = 0.5603 corresponds to A = 2.0°. 100 Monte- 
Carlo simulations are used to obtain the mean SINR and SOI power estimates. It is worth noting 
that both RCB and DCRCB allow JVto be less than the number of array elements M. We use N = 
6 and N = 8 is used for the N < M case in this example. For DCRCB, when the condition 

Re(a"a) <M-e/2 

is not satisfied, the SOI power estimate is calculated by 0 • Note from Figure 7 

that the convergence properties of both methods are quite good and somewhat similar. Since the 

errors between «■ and R can be viewed as equivalent steering vector errors, e should be chosen 
larger than eo, especially for small N. 
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[0061 ] The performances of RCB and DCRCB is now compared when the power of the 
fourth signal is varied. As in the previous example, e 0 = 0.5603 corresponds to A = 2.0°. The 
INR in Figure 8 refers to the ratio between the 4th signal power and the noise power. Note from 
Figure 8(a) that the SINR of DCRCB is much better than that of RCB when e = 0.6. However, 
when e is large, for example when e = 2.0 as in Figure 8(c), and when the INR is comparable to 
the SNR of the SOI, DCRCB has lower SINR than RCB. From the diagonal loading levels of the 
methods shown in Figures 8(e) and 8(f), it is interesting to note that the diagonal loading level of 
RCB when e = 2.0 is about the same as that of DCRCB when e = 0.6. As a result, the SINR and 
SOI power estimate of RCB when e = 2.0 are about the same as those of DCRCB when e = 0.6. 
Note also from Figure 8 that when the INR becomes close to the SNR, there is a performance 
drop in the array output SINR. One possible explanation is that when the INR is much smaller 
than the SNR, its impact on the SOI is small. As the INR increases, it causes the SINR to drop. 
As the INR becomes much stronger than the SNR, the adaptive beamformers start to form deep 
and accurate nulls on the interference and as a result, the SINR improves again and then becomes 
stable. 

[0062] Next, the case where the SOI power varies is considered, e = 0.6 is chosen and three 
cases are considered: A = -1.0°, -2.0°, and 3.0°, with the corresponding e 0 being 0.1474, 0.5939, 
and 1.2289, respectively. Figure 9 shows that as long as e is greater than Co and the SOI SNR is 
medium or high, the SOI power estimates of RCB and DCRCB are excellent. Their SINR curves 
are also quite high, but they drop when the SOI SNR approaches the INR of one of the 
interfering signals. One possible explanation is given at the end of the previous paragraph. 
Figures 9(e) and 9(f) show that when e is smaller than eo.the performances of both RCB and 
DCRCB drop drastically as the SOI SNR increases from moderate to high. This is because the 
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SOI is suppressed as an interference for this case. Note that RCB and the DCRCB according to 
the invention significantly outperform the SCB in ShNR and SOI power estimates. 

[0063] Finally, in the last two examples 100 Monte-Carlo simulations are used to compare 
the statistical performances of RCB and DCRCB. First, the true arrival angle of the fifth signal is 
uniformly distributed between 38° and 42° but it is assumed to be 40°. Figure 10 compares the 
SINR and the SOI power estimates of RCB and DCRCB obtained in the 100 Monte-Carlo trials. 
Note that the SINR mean of DCRCB is about the same as that of RCB but the SINR variance of 
DCRCB is much smaller than that of RCB (especially when e = 0.6, which is quite tight since 0 
<e 0 ^ 0.5939). Hence this example shows again that DCRCB may be preferred over RCB when 
higher array output SINR for waveform estimation is needed. On the other hand, the bias and 
the variance of the SOI power estimates of RCB are smaller than those of DCRCB. This is 
especially so for large e; note that a large e is not a problem here since the interfering signals are 
quite far away from the SOI. Hence this example also shows that RCB may be preferred over 
DCRCB in applications requiring accurate SOI power estimation including radar, acoustic, and 
ultrasound imaging. 

(0064] Next consider an example of array calibration error that consists of perturbing each 
element of the steering vector for each incident signal with a zero-mean circularly symmetric 
complex Gaussian random variable with a variance equal to 0.1. The perturbing Gaussian 
random variables are independent of each other. The calibration error is not scaled or normalized 
in any way and hence 0 <e<> <°°. In Figure 1 1, the means and variances of the SINR and SOI 
power estimates of RCB and DCRCB are compared, as functions of e. The figure shows once 
again that with a reasonable choice of e, the DCRCB according to the invention may be preferred 
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for applications requiring high SINR whereas RCB may be favored for applications demanding 
accurate SOI power estimation. 

[0065] It is to be understood that while the invention has been described in conjunction with 
the preferred specific embodiments thereof, that the foregoing description as well as the 
examples which follow are intended to illustrate and not limit the scope of the invention. Other 
aspects, advantages and modifications within the scope of the invention will be apparent to those 
skilled in the art to which the invention pertains. 
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